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Abstract. We analyze the time evolution of simple nuclear rotational wave packets 
(WP) called circular, linear or elliptic, depending on squeezing parameter rj, assuming 
that E = huiQ 1(1 + 1). The scenario of fractional revivals found by Averbukh and Perel- 
man is adapted to symmetric WP and compared to that which holds for asymmetric 
WP. In both cases various shapes are identified under these lines in particular many 
cases of cloning. 'Mutants' WP are found most often. Finally the time evolution of a 
WP formed by Coulomb excitation on ^^^U and calculated by semiclassical theory is 
also presented. 

PACS: 03.65.Sq, 03.65. Ge, 21.10Re , 21.10Ev 



1 Introduction 

It is somewhat paradoxical after many years of confirmation of quantum me- 
chanics that new features can be learned on two of the simplest systems that are 
taught in the elementary courses: the hydrogen atom and the rigid rotor. It was 
indeed realized only in the early 90's j^, ^ that a wave packet in the hydrogen 
atom evolves in time according to a universal scenario discovered theoretically 
by Averbukh and Perelman . For short times a well built wave packet moves 
around a Kepler ellipsis, and then it spreads along the ellipsis. The surprise was 
that the wave packet is seen at specific times as a superposition of fractional 
wave packets concentrated around specific points on this trajectory. This be- 
haviour results from the unitary quantum evolution and from self interferences. 
By changing the initial wave packet one can also obtain a radial wave packet 
without an analogue classical trajectory which will also experience the scenario 
of fractional revivals [Q . An extension of these results has later been been made 
by Bluhm and Kostelecky who showed the existence of superrevivals for longer 
times. 

The aim of this article is to describe the similar steps for the rotation of 
symmetric rigid rotors with possible application to nuclei and to molecules. The 
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theory of ref. |^ is exact for all times if the energy spectrum is quadratic in one 
quantum like the infinite square well Q or in two quantum numbers like |^] or 
the rigid top with axial symmetry. In ref. we studied various sets of angular 
momentum coherent states which represent an heteronuclear diatomic molecule 
or a rigid body with axial symmetry. Our purpose is to show the rotation of a 
symmetric object described by a quantum pure state in the most obvious scheme 
(i.e. if the law E = {h^ /2J){I{I+1)) holds exactly) for a wave packet specifically 
designed to rotate and what the motion is, if the initial preparation of the wave 
packet is different. Because of the quantum spreading a quantum WP does not 
rotate at all like a classical rigid body, its spreading differs according to the initial 
conditions of its preparation. The introduction of this spreading is necessary if 
one introduces the concept of a WP which is not a common practice in the field of 
rotational nuclei or molecules. We must however acknowledge that this concept 
was first introduced in the past in the reference book by Broglia and Winther Q 
on heavy ions collisions, then in a more extensive work by Fonda et. al. ||l0|| who 
have discussed the mechanism of generation of several sets of rotational coherent 
WP. Our work completes these two works in the sense that we give a full study 
of the time evolution of the WP which is now possible after the work of . Our 
WP depend also on a new squeezing parameter ij introduced in which enables 
to change the geometry of the WP in a most significant manner. 

Before proceding to our work we must also quote most elegant papers by 
Berry and Golberg and Berry who have studied the full time evolution 
of the propagator of a nuclear spin with a hamiltonian of the form L^/(2J). 
More recently various Talbot effects were also discovered with the help of renor- 
malization techniques . The occurrence of fractal dimensions both in space 
and time during the evolution of uniform WP in boxes of arbitrary dimensions 
has been also discovered . These last results show that the spreading of WP 
in a bound system is extremely rich. Our purpose is to emphasize this richness 
for two dimensional WP on the sphere. 



2 Asymmetric and symmetric squeezed WP 
2.1 Asymmetric WP 

In ref. Q we have studied the properties and the time evolution of the following 
WP called the exponential coherent WP which depends on two parameters N 
and rj. This normalized state is defined as 

* (61, 0) = J e^ '^'^ "(^^ '^+''' '^'^ '^K (1) 

y 47rsinh(2iV) ^ ^ 

For large enough values of N the probability density associated to it is symmetric 
around Ox and is peaked within a solid angle !7 — 47r/(4A^ + !)• Thus, large 
values of N correspond to defining the orientation of the axis of symmetry of 
the system in a very sharp way, quite independently on the value of the average 
angular momentum L. For the components of L one can see that 



(i.) = {Ly) = , 



(2) 
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{L,)^V{N~^). (3) 

Equation (||) is correct only for large N, for small N there is a small correction 
term not important for our discussion. Therefore, the average of L is fixed by 
the product r]N. Finally we have proved |^ that the state (|l|) is made of a 
special class of angular momentum states called in the literature by the name of 
intelligent spin states (for real values of 77) . They obey the equation 

{L^ + i'nLy)<F„+ ^ . (4) 



From this equation one can derive that the parameter 77 has the meaning of a 
squeezing parameter since 

hp = 4i (5) 



For all real values of 77 the minimum uncertainty condition is satisfied 

ALlALl = \{L^f . (6) 

The states (^ are such the analogous of the squeezed states of the electromag- 
netic field. In this article we will distinguish two limiting values of ry: the case 
77 = corresponds to a purely real WP initially without momentum which we 
call the linear WP, the case ?7 = 1, called a circular WP, is an eigenstate of L+. 
The other values correspond in the same terminology to elliptic WP. Indeed 
these WP can be derived from Gaussian WP which are respectively associated to 
linear, circular and elliptic trajectories if they evolve in 3D harmonic oscillator. 

The case of complex values of 77 does not bring any interesting freedom ||l7|| , 
the main difference lying in the fact that the equality (||) is violated. Therefore 
77 will take only real values in the following. 

The spherical harmonic expansion of the state (1) is written in terms of 
coefficients 5/m as 

S'.-E^^a^^mI^^,^), (7) 

IM 



2N 



'JIM = 



X 



sinh(27V) 



w 



(8) 

(1 + 7/)'(l - 77)'' {ll'00\IQ){Wl - l'\IM) 



For 77 = 1 only the 6's with M = I are non-zero. This is why these wave packets 
are called circular. In the case 77 = it is simpler to make a rotation of the 
coordinate axis to get an eigenstate of with eigenvalue 0, the expansion of tf^o 
takes the form 

^^ = Y.bjoY^{e,c^) . (9) 
with 6/0 given in terms of spherical Bessel function of the first kind by 



&/0 = 



y(27TT)y^/,+i/2(iv). (10) 



sinh(27V) 
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2.2 Symmetric WP 



The situation we want to describe is that of an even, deformed nucleus with an 
axis of symmetry such that (Lz) — K — and with a plane of symmetry 
perpendicular to this axis. No other degree of freedom will be included. We 
therefore consider instead of (|l|) the following symmetrical combination 



9± 

2 



^A^ sin 0{cos 4^-\-iri sin 0) 



-N sin 6{cos 0+i?7 sin 0) 



(11) 



Clearly for large N it represents two antipodal waves, within the same solid angle 
as above and with the same symmetry around Ox. The normalization constant 
C+ depends now both on TV and rjN because of an interference term between 
the two waves. It is given by 



2n 



sinh 27V sin 2riN 



2N 



(-1/2) 



(12) 



It should be stressed that the main difference in the time evolution of (|l|) and 
( pr] ) comes from the interference between these two parts. As we will see it will 
become much different from a normalization factor. 

It is straightforward to verify that eqs. (^) to (^ are still valid. We now 
express the wave packet in terms of new coefficients 6/A/+ as 



1^, 



E 

IM 



nM+ YMiS, 



(13) 



and these coefficients are related to the bjM by 



5/M+ 




3/M 



V2b 



IM 



for / odd 
for / even. 



(14) 



3 Time evolution 



3.1 The rigid rotor assumption 

Let us study the evolution of if",,-!- for rj — and 77 = 1 assuming that the energy 
levels obey the rule 

/?, = ^/(/ + l). (15) 

In the following we will use t^o = fi/(2J) and we define an average angular 
momentum by 

/(/+l) = (*,+ |i>„+). (16) 

We can now use the work of |^ to analyze the structure of the wave packets at 
time t knowing simply that the wave packets are given by 

1r,+ (0, 0, i) = ^ bjM+ Ylj{e, (p) e--° ^(^+1)* (17) 

IM 
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It is necessary to introduce two characteristic times as in 

Tci = — 77TTT' -^'■cv ^ — ■ (18) 

oJo 21+1 ujQ 

Note that T^v >> Tci if I is large enough. We have with Q and (|l|) 

^„+(0,0;Trev)=>Z'r,+ (0,0;O). (19) 

3.2 The outline of the theory of fractional revivals 

Let us sketch the general theory of |^ which allows to find a universal scenario 
of fractional revivals. Near t = (m/n)Ticv, where m/n is an irreducible ratio of 
two integers, it has been proved in |^ that 

^ ^ ^ a,, (20) 

s=0 

This sum contains 9 = § non-zero terms if n is even and q = n terms if n is 
odd. The integer I = j ii n is a. multiple of 4, Z = n in the other cases (therefore 
the difference between I and q is necessary to take into account that the Ug with 
even s are zero if n is even and not multiple of 4). The modulus of Ug is simply 
for all the cases where it is non-zero equal to: 

kl = ^ (21) 

The phase of as a function of m, n and s can be found in an older paper 
written in a quite different context . Therefore one writes after the insertion 
of into (0) taken at time t = f T^ev 



/-I 



n ^ — ^ 



s=0 
i-1 



.IM 



(22) 



- Y.as^^^\e,<p-U) . (23) 

s=0 



is) 

Any wave packet is a sum of q fractional wave packets if'^j , each with a different 
effective time ts 

Tll> S 

= (- + y)Trev. (24) 
The fractional wave packet is defined as 

W^l^ {0, <P;t,)^Y. Yi, e--"^*^ . (25) 

IM 
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3.3 Clones and mutants 

There are several cases for which tf'ci is a clone of the initial ^ defined by or 

i) if the expansion contains only biM with the single values M = I, this is 

the case for the wave packet (|l|) for 77 = 1. One obtains for each tg 

^ci{e,(f>;ts) = i'„{e,4>-^ots) ■ (26) 

This result is independent of the distribution of the 6//'s, it applies to all 
kinds of circular wave packets. The sum (|2^) is a superposition of q clones 
of the initial wave packet. In |6j this conclusion was derived for an infinite 
one-dimensional well. In our case with rj = 1 the q clones are localised 
around q axis forming a symmetric fan in the Oxy plane. 

ii) if LOots is a multiple of 2tt, i.e. if is non-zero for the value 

So = n — m , (27) 

then 

'Z^ci(0,<^;tso)-*,(e,0,O) . (28) 

This occurs when n is odd or even and not multiple of 4. For these times 
there is a clone at the place of the initial wave packet with the amplitude 
ttsg ■ The existence of this clone is independent of the 6/m and the parameter 
V- 

In 1^ we have proposed to call mutants the fractional wave packets tf'ci which 
are not clones and which have a different geometry due to the phase e~'"°^*= in 
(p5|). We have also described the change of the shape of mutants as a function 
of r] for fixed ts ■ This is made by considering two extreme cases where the shape 
of tf'ci is simply understood. For 77 = 1 tf/ci is a clone of the initial wave packet 
located around a given direction in the Oxy plane. For ?/ = the expansion 
shows that wave packet !Po a-t ^ = as well as for all t and also each 'F^ 
have cylindrical symmetry around Ox. S'ci is therefore peaked around one or two 
rings around Ox. For intermediate values of rj the mutants have shapes making 
a smooth transition between the two extreme cases. 

This theory fits perfectly to the case of the asymmetric WP which contains 
terms with even and odd values of / in equations @) , (^0|) and (^2|) . In the case 
^ = ^ eq. gives tg — T^-cv, Os = Sgi and the WP is periodic with period 
ts = rrcv/2. This reduced periodicity can be seen immediately from the fact 
that /(/ + 1) is always even. It is therefore enough to consider times such that 

m ^ 1 
n — 2- 



3.4 Additional rules for symmetric WP 



In the symmetric case corresponding to (11) and ( [13| ) with (|14|), an additional 
symmetry arises from the fact that eq. (20) is now used for even values of /. If 
we put n = 4n' in (Effl we obtain q = I = 2n' values of s in the general case. 
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However if / is always even I takes the value of n'/2 if n' is a multiple of 4 or 
n' is odd or even and not multiple of 4. For n' even and not multiple of 4 the 
Qs with even s are zero. Therefore the number q of fractional waves is largely 
reduced and becomes 

_ ( n' for odd n' , , 

1 n'/2 for even n' . ^ ' 



For n' = 1 one can write 



for n' = 2 



g-2i.7(/+l)/4^gi.//2 ^ (30) 



g-2i7r/(/+l)/8 ^ g-i37r7/4 _ ^gj^^ 

Inserting these values in ( [171 ) for the particular case of circular WP (77 = 1) one 
obtains the following periodic behaviour 

tf'i+(0, dp, r,.ev/4) = ^i+{e, + ^/2, 0) , (32) 

<f'i+(0,0,T,ev/8) = 'f'i+(f?,0-37r/4,O) . (33) 

For 77 ^ 1 the insertion of the phases (^9|) into (|l^) produces a single WP different 
from the initial one, i.e. a single mutant. For rj = 1 the time interval necessary 
to study the evolution is t = Tj-ev/S, but for ry 7^ 1 it is still t = T].ev/2. 

At fractional times of Ticv/4, like t = Tjcv ^ the following formula holds for 
r] — I and describes fully the cloning of these species of symmetric WP 

I 

^MO, 0, r,.ev £7) = E - T(y + 1^), 0) , (34) 

with the general result that the sum contains q — n' /2 non-zero terms for even 
n' (with I = n' /2 if ?i' is a multiple of -i, I — n' in other cases and as is zero if n' 
is even and not multiple of 4 and s is even) . 

For example the application of this formula for t — Trcv/lS gives two clones 
oriented along perpendicular directions in the Oxy plane since 

tf'i+(0,<^,T,ev/16) = 00*^1(0, <^-7r/8,0) + 01*^1(0,0-5^8,0) . (35) 

Similarly, for t — T^cv/2A one obtains 3 clones etc... For 77 7^ 1 the clones should 
be replaced by a corresponding number of mutants. 

If on the other hand n is even and not multiple of 4, i.e. if n = 2?t,' with 
n' odd, the limitation of even / in eq. (20) leads to I = n' therefore all the 



are non-zero but should be calculated with an effective ratio 2m /n' instead of 
m/n. The number of fractional WP is then the same as for the asymmetric case. 
Similarly, if n is odd there are q = n fractional WP. Note that for odd n and 
even n not multiple of 4 eq. (|7|) and (|8|) also hold: one of the fractional WP is 
a clone at the place of the initial WP. 
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4 Numerical calculations for the exponential coherent 
WP. 

4.1 Comparisons between the symmetric and asymmetric WP 

As already stated above, the WP ( pd] ) differs from (|l]) by the interference terms 
which are almost negligible if one calculates the probability density at t=0. 
However the previous discussion on the time evolution shows a large number of 
differences. Therefore a comparison between asymmetric and symmetric WP is 
presented in the figures of this letter. With high enough TV one obtains a generic 
behaviour. For ry = 1 we have chosen TV = 14 i.e. / = 13.5. This value of I 
corresponds to the average excitation energy of the K ^ {) rotational band of 
^•^^U that is obtained from Coulomb excitation by 200 MeV ''"Ar projectiles. 
Coulomb excitation does not produce weights hiM so simply distributed to give 
a WP like 

Fig. |l| shows 'carpets' for the two cases with TV = 14 and rj = 1. The carpet 
is the intersection of the plane Oxy with the probability density as a function 
of time for < i < T^^^jl. For t < Trov/SS = Tci the WP spread and occupy 
all values of phi. For t > Td there are clear windows near values {m/rijTj-ev for 
which the WP is formed of q clones, with q given as discussed in the text for 
asymmetric and symmetric case respectively. In the former case one has large 
windows for m/n = 1/4 {q = 2), m/n = 1/6 and 1/3 {q = 3), m/n = 1/8 and 
3/8 {q = 4) and it is possible to see as much as g = 8 clones without much 
interference between them (see Fig. El). In the latter case we obtain instead due 



to eqs. ( |3l| ) and (33) a slide angle of — 37r/4 since the following equation takes 
place 

(0, 0, t + r,,v/8) = If 1+ (0, - 3^4, t) . (36) 

Therefore the interval Trov/S is separated into four stripes into which essentially 
only a single large time window is seen for m/n = 1/16 with q = n' /2 = 2 clones. 
By zooming within Ti.ov/8 one can see other windows as well as we will see in 

Fig- 1 _ _ 

Fig. H and y show the probability density as a function of 9 and (j) for a 
selection of times where one can underline similarities or differences between 
the two cases. For < i < Tci the WP essentially spreads in cj) because the 
corresponding values of n produce too many clones with big interferences 
For t = Trov/50 the two peaks forming the symmetric WP interfere strongly 
and the WP occupies already all the values of (p. In the next sequences of times 
one identifies {q,q') clones respectively in the asymmetric and symmetric cases: 
m/n = 1/16 (8,2), m/n = 1/3 (3,3), m/n = 1/24 (12,3): there the 12 clones 
are strongly interfering, m/n — 1/5 (5,5), m/n — 1/6 (3,3). 

Fig. 1^ shows the carpets corresponding to r/ = with TV = 110 and not 
TV = 14. The distribution of the bjM with / depends indeed very much on the 
value of rj. Using TV = 14 with r] ~ produces too low value of / and a rather 
poor carpet. The value TV = 110 corresponds to / = 10 and leads to richer 
figures. For this value of rj the WP is axially symmetric for all times, the carpet 
is therefore defined as the quantity 2tt s'mO \ipo{9,t)\'^ as a function of 6 and t. 
Both carpets are also symmetric around t — Tj-cv/i. However only the carpet 



Pictures of quantum nuclear rotation beyond the correspondence principle 



9 



showing the symmetric case has an additional symmetry around 9 = Tr/2. Indeed 
a time r which fulfills for all I the equation 

exphic^or/(/+l)] = (-1/ (37) 

does not exist. The maxima of the Wj^^\0,ts) occur for definite values of 9 and 
are not always distributed symmetrically around 7r/2. These fractional waves 
can be understood as waves on a sphere with maximum along a ring j^. For 
t = Trev/4 there is a unique large ring at 9 = n/2. Moreover, since 6/o and Yjq 
are both real, values s and s' may exist such that 

'F^j'^^'F^f*, (38) 

for example for t — 1/10 Tiov there are 5 fractional WP with < s < 4. For 
So = n — m — A there is a clone such that ^j.f'^ is equal to the initial ^. This 
WP is represented as a dark dot at 6* = in the left part of Fig. ^ and there are 
two pairs of fractional WP, one for s = and s' — 4, the other for s = 2 and 
s' = 3. By addition of these 4 WP one obtains only two rings represented as two 
other large dots for different 9. In the symmetric case this number is doubled for 
Tj-ov/lO. As discussed above the existence of the clone relies upon valus of m/n. 

In the right part of Fig. ^one can most often detect a ring a,t 9 — Tr/2. Similar 
carpets have been discussed in [p6| for the infinite well in one dimension. 



4.2 Geometrical shapes of clones and mutants 

We present in Fig. || a short selection of shapes of WP for 4 different values of 
r] and the same time t — Trev/8. Only the part of the jlZ'p with 2 > is shown, 
note that (j) is defined as the angle with the Oy axis. In the asymmetric case one 
obtains 4 clones for 77 = 1. When 77 is decreased these clones become 'mutants' 
which develop a larger spread in the angle 9. For 77 = 0, as well as for rj = 1/25, 
one obtains two rings around the symmetry axis Oy as well as a maximum along 
this axis. In the symmetric case we illustrate for the circular WP the law given 
by eqs. ( ^ , ^ : there is a periodicity with a shde angle — 37r/4. For different rj 
the clone is transformed into a single mutant with a larger spread in 9. For the 
case of the linear WP this mutant is cylindrically symmetric around Oy. Note 
that the mutant has, for 77 = 1/25, an interesting dissymmetrical shape which is 
necessary in order to make the transition between the clone for rj = I and the 
two rings for 77 = 0. 

The WP studied in this letter exibits well defined shapes and a well under- 
stood evolution. This is due to the assumptions that the law (|l^) holds exactly 
and that the 6/ m are simply distributed according to (||) in order to correspond 
to the simple exponentials (|l|) or (0). 

5 Calculations with realistic WP 

It is quite possible to improve our study in order to incorporate ingredients which 
make it more realistic. It is well known H |l3 that during Coulomb excitation 
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(CE) a deformed nucleus is excited to a coherent mixture of rotational states. 
This superposition is also peaked around a mean value of angular momentum, 
so one can expect similar features as predicted by scenario of Averbukh and 
Perelman. The most clear case is CE with backscattering as in this case excited 
WP has cylindrical symmetry (only Yq components, i.e. 77 = or linear WP). 
The partial waves and fractional revivals have then topology of rings on a sphere. 
Therefore for presentation of shapes only one angular variable {9) is sufficient. In 
Fig. ^ we present the time evolution of WP obtained by CE of ^ssy bombarded by 
^°Ar at E=170 MeV. The amplitudes of excitation of given / angular momentum 
eigenstates have been calculated within semiclassical theory of CE The left 
part shows the 'ideal case', i.e. when Ej's follow perfect rotor dependence /(/+1), 
the right corresponds to time evolution with energies taken from experiment. In 
the latter case it is necessary to define time scales of the evolution. In general 

^-^wk~i' ^''^wtri- ^^^^ 

For perfect rotor case these times are determined unambigously: Tjcv ~ (27r/?i^) J, 
Tf-i = Ti.ev/(2/+ 1). For experimental energies we adopt the following estima- 
tion for Trev The energies of the ground state rotational band for ^ssy ^^^.g 
approximated by the two parameter formula 

£;/~a/(/ + l) + &[/(/+l)]2 (40) 
with a = 7.5 MeV, b = -0.004 MeV pi. Our 'estimated' T^v, used for the 



figures (there is no exact revival at all), is then defined by (39) and ( |40| ) taken 
for / determined from the equation 

/(/+1) = (L2) (41) 

i.e determined by coefficients of CE rotational WP (see Fig. 8, p. 82 of rcf. [||). In 
Fig. 1^ we display shapes of WP (precisely the probability density as a function of 
9 angle) for short times, showing spreading of the WP, initially localized at poles, 
and clear few rings structure of WP at particular times (3/400 and 4/400*Ti.cv)- 
In Fig. ^ we show the same WP after much longer evolution, exhibiting structures 
very similar to those observed for very short times. 

Although the 'carpet' for experimental energies isn't as regular as that for 
'ideal' ones, still strong revivals of WP occur. The absolute time scales for nu- 
clear rotation are very short (Trcv ^ 10~^^s, Td ^ 10^^°s) and are still beyond 
time resolution of present experimental techniques. It is clear that the picture 
which emerges from this study underlines that the time evolution of nuclear WP 
depends strongly on their excitation mechanism as well as from the assumption 
of a rigid body value for the moment of inertia. The inclusion of more realistic 
ingredients introduces less symmetry and less regularity in the time evolution. 



6 Conclusions 



There is no doubt that the time evolutions described in our article have been 
calculated by observing strictly the rules of quantum mechanics for a pure state. 
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It is reasonable to believe from the works of |^ and |jT3 that a WP formed by 
Coulomb excitation should move as described in section l5|. However there is not 
much hope that the tools of experimental nuclear physics can evolve in such 
a way that this time evolution is observable since one needs to probe the WP 
after its formation with the help of a second projectile. In future the creation 
of nuclear rotational wave packets and observation of their evolution should be 
in principle possible with gamma lasers. Nevertheless the picture of the rotation 
one should have in mind must be that given above, i.e. a picture where the 
interference effects play a dominant role and differ according to the symmetry 
of the rotor and its initial preparation. 

The field of molecular physics looks more appropriate for observation and 
many theoretical calculations as well as several experiments have already been 
realized. Rotational coherence effects of large molecules have been studied the- 
oretically and observed for about ten years [|9|-||2^ but are generally under- 
stood as thermal average of rotational quantum beats. The rotation of a single 
molecule of O2 |Q or HB-DC |23| within a supra molecular bearing has been 
recently observed. The rotation is induced thermally and there is no reason to 
believe that the system is in a pure state. In the frame of optimal control theory 
[ p4| one has also shown recently how to create rotational coherent WP of 
asymmetric-top molecules by designing tailored microwave pulses |2^, ^ . The 
order of magnitude for the time scales of molecular rotations varies from 10~^^s 
[ p7| to 10~^^s 1^ and is available for experiment. However, the known excita- 
tion mechanisms lead mainly to rovibrational states and an excited superposition 
is usually composed of smaller number of components than rotational coherent 
states discussed in this paper. A most important challenge for our purpose would 
be to find the excitation mechanism that enables the generation of our squeezed 
WP. Excitation mechanism of a rigid rotor with intense laser field is still under 
study but like Coulomb excitation of nuclei, this mechanism leads only to linear 
waves (with ?/ = 0). 

The study of the time evolution of WP in simple sytems produces behaviour 
which is full of surprises. A recent work on WP transmitted through a potential 
barrier has detected the presence of high momentum components which are 
present during a short time |^ an event which is also a challenge for detection. 
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Fig. 1. Time evolution of circular, asymmetric (left) and symmetric (right) WP 
with iV = 14 in the 'carpet' representation. The |'f'(6', 0, i)p is presented in the 
contour plot for fixed 9 = 7r/2. 



Fig. 2. Shapes of circular, asymmetric WP in the initial stages of the evolution 
and some fractional revival times. Note different vertical scales used in order to 
show details. 



Fig. 3. The same as in Fig. but for symmetric WP. 



Fig. 4. 'Carpet' representation of the time evolution of linear, asymmetric (left) 
and symmetric (right) WP with iV = 110 and 77 = 0. In this case the probability 
density 27rsin6'|if'(6',t)p is presented in the contour plot. 



Fig. 5. Shapes of clones and mutants for asymmetric (left) and symmetric (right) 
WP with iV = 14 and different -q at time t = Trov/S. The |!f (6*, (j), t)p is presented 
in spherical coordinates as the radial coordinate, 9 is an angle with respect to 
Ox axis which is the axis of cylindrical symmetry in this case. (Only upper part 
of the probability density is displayed) 



Fig. 6. Time evolution of nulear rotational wave packet obtained in CE of ^ssy 
presented in 'quantum carpet' representation. Contours of 2-Ksm.9\]/\^ are plot- 
ted. 



Fig. 7. Shapes of nulear rotational wave packet obtained in CE of ^ssy during 
a short time evolution. 



Fig. 8. The same as in Fig. ^ for some particular much longer times. 



N=14 r|=l all I N=14 r|=l even I 




Fig. 1. Time evdlution (if rirculiu". aaymmetric (left) and aymmetric (right) WP 
with = 14 in the 'rarpet' repreaenttitidn. The {^(0: i^: *)|^ ia presented in the 
nintdur plot for fixed 8 = 7r/2. 
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Fig. 2. Shapes nf rimilar, asynunietrir WP in the initial stages nf the evnhitinn 
and snme firartinnal revival times. !Nnte diffifTHit vertiral scales iLsed in nnder to 
shnw details. 
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N=14 T|=l even I 




Fig. 3. The sanie as in Kg. 2 biit fihr symmetrir WP. 



N=110r|=0 all I N= 1 1 r|=0 even I 
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Fig. 4. 'Carpet' repreatHittitidn (if the time evdlution (if linear, aaymmetric (Itift) 
and aymmetric (right) WP with ^ = 110 and 17 = 0. In this rjiae the pniba-hility 
dtinfdty 27rsiinfl|!P'(fl. *)|^ ia preatinttid in the ("j(int(iur pldt. 



N = 14 t=1/8*T 



all I 
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Fig. 5. Shiiptai of rltrntsi and mutimtH for aaymmetric (left) and aymnmtrir (right) 
WP with ^ = 14 and different 17 at time * = T„^/B. The \W(8,^,t)\^ is presented 
in apheriral cjotirdinatea aa the radial cjinH-dinate. ^ ia an ang^e with reapert to 
Ox axis whidi ia the axia (rf cylindriral aymmetry in thia raae. (Only upptr part 
(if the pnibahiUty denaity ia diaplayed) 



E|=l(l-i-1) Exp. energies 
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Fig. fi. TiiDB (^fnhiticin nf mileax mtalinnal waifE pajcliEt obtained in CE nf 
presented in 'quantum carpet' representation. Ccmtoiirs nf 25rsin^|!P|^ are plot- 
ted. 




Fig. 7. Shapes of nulear rotational wave packet obtained in CE of ^^^U during 
a short time evolution. 
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Fig. 8. The same as in Fig. 7 for some particular much longer times. 



